Interaction and disorder in bilayer counterflow transport at filling factor one 
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We study high mobility, interacting GaAs bilayer hole systems exhibiting counterflow superfluid 
transport at total filling factor v = 1. As the density of the two layers is reduced, making the bilayer 
more interacting, the counterflow Hall resistivity (p xy ) decreases at a given temperature, while the 
counterflow longitudinal resistivity (p xx ), which is much larger than p xy , hardly depends on density. 
On the other hand, a small imbalance in the layer densities can result in significant changes in p xx 
at v = 1, while p xy remains vanishingly small. Our data suggest that the finite p xx at v = 1 is a 
result of mobile vortices in the superfluid created by the ubiquitous disorder in this system. 

PACS numbers: 73.43.-f, 71.35.-y, 73.22.Gk 
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Interacting bilayer systems of two-dimensional carriers 
in the limit of zero tunnelling can exhibit superfluidity 

in a peculiar "counterflow" transport configuration, 
where currents of equal magnitude are passed in oppo- 
site directions in the two layers. This phenomenon occurs 
when the bilayer is subject to a perpendicular magnetic 
field (B) so that the total Landau level filling factor, v, is 

1 (layer filling factor 1/2) and the inter-layer interaction 
is sufficiently strong to stabilize a quantum Hall state 
(QHS) Q]. The physics of this QHS can be understood 
by pairing the particles in the lowest Landau level of one 
layer with the vacancies in the other layer, hence forming 
neutral objects (excitons) which condense at the lowest 
temperatures Q- A similar superfluid was predicted to 
occur in a closely spaced electron-hole bilayer at B = 
0. Key experimental evidence for the particle- vacancy 
pair formation at v = 1, and the ensuing condensation, 
has been provided so far by two types of experiments. 
Inter-layer tunnelling measurements have shown a much 
enhanced tunnelling conductivity in the v = 1 QHS 

a behavior reminiscent of a Josephson junction's. Most 
recently, counterflow transport measurements have re- 
vealed that both the longitudinal (p xx ) and Hall (p xy ) 
resistivities vanish in the v = 1 QHS in the limit of zero 
temperature (T — > 0) ■ The vanishing of p xy is espe- 
cially important since it directly demonstrates that the 
counterflow current is carried by neutral particles, that is 
particle-vacancy pairs which have zero electrical charge 
and therefore experience no Lorentz force. 

An outstanding puzzle is what causes the dissipation 
in the currently available samples and why experimen- 
tally there is no finite critical temperature below which 
the counterflow dissipation vanishes as the theory pre- 
dicts 0]. Here we study the counterflow transport in 
strongly interacting, high mobility GaAs hole bilayers as 
a function of total density (p) as well as layer density im- 
balance. Our data show that, at a given temperature, the 
v = 1 QHS counterflow p xy decreases when p is reduced 
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to increase the inter-layer interaction. But the counter- 
flow p xx , which is always much larger than p xy , barely 
depends on density. Furthermore, small changes in the 
layer densities can substantially change the counterflow 
p XXl while p xy remains vanishingly small. The vanishing 
p xy at low temperatures demonstrates that inter-layer 
interaction is responsible for the particle-vacancy pair- 
ing. The counterflow dissipation, signaled by the finite 
p xx at finite temperature, has been attributed to the ex- 
istence of disorder-induced mobile vortices which move 
across the superfluid current, much like in superfluid he- 
lium 0, H, @ • Our observed dependences of counterflow 
p xx on total density and layer density imbalance are con- 
sistent with this picture. 

Our sample is a Si-modulation-doped GaAs double- 
layer hole system grown on a GaAs (311) A substrate. It 
consists of two, 15nm wide, GaAs quantum wells sepa- 
rated by an 8nm wide AlAs barrier. The top and bottom 
barriers are Al0.2Ga0.sAs layers. We used a Hall bar ge- 
ometry of 100/im width, aligned along the [Oil] crystal 
direction. The Hall bar mesa has two current leads at 
each end, and three leads for measuring the longitudinal 
and Hall voltages across the bar. Diffused InZn Ohmic 
contacts are placed at the end of each lead. We use a com- 
bination of front and back gates to selectively deplete one 
of the layers near each contact, in order to realize inde- 
pendent contacts to each layer [l(j. As grown, the densi- 
ties were pr = 3.1 x 10 10 cm~ 2 and ps = 3.8 x 10 10 cm~ 2 
for the top and bottom layers, respectively. The mobility 
along [Oil] at these densities is approximately 32 m 2 /Vs 
[rH . Top and bottom gates were added on the active area 
to control the layer densities. The measurements were 
performed down to T = 30mK, and using low-current 
(0.5nA-lnA), low-frequency lock-in techniques. 

In our counterflow measurements, two leads contacting 
opposite layers at one end of the Hall bar are used to drive 
a current in and out of the sample, while the leads at the 
other end are shorted so that the same current, but in 
opposite directions, flows in both layers [Fig. 1(b) inset]. 
We define the resistances in the counterflow configuration 
as the corresponding voltage drops along or across the 
Hall bar measured in one single layer, the bottom layer 
in our case, divided by the current flowing in each layer. 
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FIG. 1: (a) p xx and (b) p xy counterflow resistivities vs 1/v for different values of d/ls, corresponding to different total bilayer 
densities, (c) Counterflow p xx (hollow symbols) and p xy (filled symbols) at v — 1 vs T -1 for different d/ls- The data of all 
panels were taken when the two layers have equal densities. Inset in (b) illustrates the geometry for counterflow measurements. 



The longitudinal resistivity p xx is the measured longitu- 
dinal resistance divided by two, the number of squares 
between the voltage probes in our sample, while the Hall 
resistivity p xy is simply the Hall resistance. These resis- 
tivities revert to those of a single layer when the coupling 
between the two layers is negligible. In our counterflow 
measurements a small fraction of the current injected in 
one layer leaks unintentionally into the opposite layer. 
The leakage current at v = 1 is about 1% of the total 
current at 30mK and increases with T to 15% at 700mK; 
this behavior suggests that the leakage is determined by 
factors (e.g. defects in the AlAs barrier) other than the 
enhanced inter-layer tunnelling at v = 1 [j| . The leakage 
translates in a slightly lower counterflow current, result- 
ing in a small error in p xx and p xy . We emphasize that 
none of our conclusions are affected by such errors. 

The phase-space of the v = 1 QHS is parametrized by 
the ratio d/ls, of the mean inter-layer spacing (d) and 
the magnetic length (Ib = yhJeB) at v = 1. This pa- 
rameter quantifies the ratio between the intra- and inter- 
layer interaction: when d/ls 3> 1 the latter is negligible 
and the physics of the bilayer is essentially that of two 
independent layers, while for d/ls — 1 the inter-layer in- 
teraction is sufficiently strong to give rise to collective bi- 
layer phenomena such as the v = 1 QHS jl2^ . Exploring 
this desired regime of small values of d/ls poses exper- 
imental challenges. Small d and/or large Ib, hence low 
densities, are needed to reduce d/ls- Reducing d, how- 
ever, increases the single particle inter-layer tunnelling, 
which impedes the fabrication of independent contacts 
and also can reduce the bilayer physics to that of a single 
layer. From this perspective, GaAs holes are desirable 
for exploring interacting bilayer phenomena: owing to 
their relatively large effective mass, the inter-layer bar- 
rier thickness can be reduced without substantially in- 
creasing the tunnelling, hence placing the hole bilayers 
in a very strong interacting regime with small tunnelling 
[l3j . Indeed the data presented here correspond to the 



smallest d/ls at v = 1 yet reported for samples with in- 
dependent contacts, also evinced by the fully developed 
v = 1 and vanishingly small counterflow p xy at relatively 
high temperatures. 

In Fig. 1 we present the counterflow p xx and p xy data 
when the two layers have equal densities (balanced). The 
different d/ls indicated in the figure (1.24 to 1.42) cor- 
respond to different p (4.5 to 6.1 x 10 10 cm -2 ). Data 
of Fig. l(a,b) show p xx and p xy vs 1/V, measured at 
T = 600mK, illustrating that p xx 3> p xy at this temper- 
ature for all values of d/ls explored here. The data of 
Fig. 1(b) further reveal that p xy at v — 1 decreases when 
djls is reduced. For example, as shown in Fig. 1(b), at 
T = 600mK p xy at v = 1 drops by more than a factor 
of 4 when d/ls is reduced from 1.42 to 1.24. This ob- 
servation demonstrates that the inter-layer interaction is 
responsible for the particle-vacancy pairing, signaled by 
the vanishing counterflow p xy 14]. Remarkably the same 
change in d/ls induces a much smaller change (~ 30%) 
in counterflow p xx at v = 1, as shown in Fig. 1(a). 

Figure 1(c) summarizes the counterflow p xx and p xy 
data at v = 1 vs T , for various values of d/ls- 
These data further demonstrate that, in contrast to p xy , 
the counterflow p xx is much less dependent on d/ls at 
any temperature. Fitting an exponential dependence 
p xy oc exp(— A#/2T) to the data of Fig. 1(c) yields a 
Hall energy gap Ah which increases from 6.8K to 9.2K, 
as d/ls is reduced from 1.42 to 1.24. The gap extracted 
by fitting p xx oc exp(— A/2T) to the p xx data is A = 2K , 
almost independent of d/ls- This finding is very reveal- 
ing: unlike the pairing, which becomes stronger when the 
bilayer is made more interacting (smaller d/ls), the dissi- 
pation is barely dependent on the inter-layer interaction 
in our sample. We will discuss this observation in our 
closing paragraph. 

In view of the large counterflow p xx and p xy anisotropy 
at v = 1, it is useful to examine the role of the mobil- 
ity anisotropy in our sample. Our measurements show 
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that the Hall drag remains close to the quantized value of 
h/e 2 as T is increased, and hardly depends on the crystal 
direction along which the current flows. The close quan- 
tization of the Hall drag at higher T is consistent with a 
small countcrflow p xy |jj and suggests that the counter- 
flow p xy does not depend on the crystal direction. On the 
other hand, our bilayer measurements on samples with- 
out independent contacts (parallel flow) show that p xx 
at v = 1 is about 2.5 times larger for current oriented 
along the low mobility direction, and exhibits activated 
T dependence with an energy gap that is independent of 
the current direction. These findings suggest that, while 
the counterflow p xx and p xy anisotropy at v — 1 may 
be partially enhanced by the mobility anisotropy in our 
samples, the latter is not the main factor behind the ob- 
served counterflow anisotropy and the strong pairing. 

Two noteworthy aspects of Fig. 1 data are (1) the large 
difference between the counterflow p xx and p xy , and (2) 
the absence of a finite temperature below which the coun- 
terflow transport is dissipationless. These can both be 
understood by considering the motion of unpaired, mo- 
bile vortices in the superfluid flow. Vortices in a super- 
fluid are subject to a Magnus force, perpendicular to the 
direction of the superfluid current, which causes them to 
move across the superflow. The unpaired vortex motion 
results in phase slip and implicitly, dissipation, hence the 
finite values of p xx while p xy remains close to 0. In the 
case of the v = 1 QHS counterflow superfluid the picture 
is slightly complicated because the vortices possess elec- 
trical charge and dipole moment 0,but the conclusions 
of the above argument remain valid Q • Theoretical stud- 
ies [H, 0| have further shown that the source of unpaired 
vortices in the v = 1 superfluid is disorder. The experi- 
mentally observed absence of a critical temperature be- 
low which the superfluid is dissipationless suggests that 
mobile vortices are present at any finite T in our sam- 
ple. As the temperature is lowered the vortices become 
pinned by disorder, form a vortex glass, and lead to dis- 
sipationless flow only in the T — limit. We emphasize, 
however, that a quantitative understanding of the role of 
disorder, vortex formation, and pinning in v = 1 coun- 
terflow superfluid is currently lacking. 

It is instructive to compare our data with similar coun- 
tcrflow measurements reported for GaAs electron bilayer 
systems 0, 0|. In Ref. @ the T dependence of p xx 
and p xy at cI/Ib = 1-48 was measured in the range of 
30-500mK and, in contrast to our GaAs hole data, little 
difference was re por ted between counterflow p xx and p xv 
in this T range [17|]. On the other hand, in Ref. jl6| . 
the counterflow p xy is reported to be much smaller than 
p xx at 40mK for a bilayer electron system with a rel- 
atively large d/ls = 1.57, qualitatively consistent with 
our hole data. While we cannot explain these similari- 
ties and differences, we mention three distinguishing fac- 
tors of holes, besides their mobility anisotropy. First, the 
larger cyclotron effective mass of GaAs holes vs electrons 
causes more Landau level mixing; whether this leads to a 
stronger pairing is not clear. Second, there is less inter- 
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FIG. 2: Counterflow p xy at v = 1 for different d/ls, vs T. 
The inset shows an empirical scaling: the p xy data of the 
main panel collapse approximately into a single curve when 
plotted vs T-B v -i, where B v= \ is the magnetic field at v — 1. 



layer tunnelling in the hole system |13J , Third, it is pos- 
sible that the holes are more fully spin polarized at v = 1 
thereby increasing the strength of the v = 1 QHS (l8j|. 

It is interesting to examine the counterflow p xy data 
at v = 1 on a linear scale and down to the lowest T, 
as shown in Fig. 2. For all oI/Ib, the counterflow p xy 
remains small up to a certain temperature at which it 
starts to increase relatively sharply. As apparent from 
Fig. 2, this onset increases as d/ls is reduced. In Fig. 
2 inset we show an empirical scaling of the p xy data at 
different values of d/ls- When the counterflow p xy is 
plotted vs T ■ B v —x, where B v —% is the magnetic field 
at v = 1, the data points collapse onto a single curve. 
Although the decrease of p xy with decreasing d/ls and 
the resulting stronger pairing is expected, the origin of 
such simple scaling is unclear. When extrapolated to 
large values of djls the scaling shown of Fig. 2 must 
break down, since it does not predict a critical d/ls for 
the disappearance of the v — 1 QHS. Experimentally 
such critical d/ls exists an d m GaAs hole bilayers 
with negligible tunnelling is between 1.6 and 1.9. 

We now turn to a study of counterflow transport at 
v = 1 as a function of layer density imbalance. Us- 
ing the top and bottom gates on the active area, we 
keep p constant while transferring charge from one layer 
to another. We define the layer density imbalance as 
Ap = (ps — Pt)/2, and measure the counterflow resis- 
tivites of the bottom layer. Depending on the sign of Ap, 
the latter is either the majority density layer (Ap > 0), or 
the minority layer (Ap < 0). In Fig. 3(a,b) we show ex- 
amples of counterflow p xx and p xy vs B traces, measured 
at T — 345mK for different values of Ap, at constant 
p = 5.45 x 10 10 cm~ 2 corresponding to d/ls = 1.35. The 
data of Fig. 3(a) clearly demonstrate that p xx at v = 1 
and therefore the dissipation, can be greatly changed by 
a small layer density imbalance: a Ap/p of ~ %5 changes 
the counterflow p xx by more than one order of magnitude. 
A similar dependence of individual layer p xx with density 
imbalance was recently reported in GaAs electron bilay- 
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FIG. 3: (a,b) Counterflow p xx and p xy vs B traces, measured at a constant total density p = 5.45 x 10 10 cm -2 (cL/Ib = 1.35) 
but for different layer charge imbalance configurations, Ap. (c) Counterflow p xx and p xy vs T -1 at different values of Ap, 
indicated by the different symbol shapes. Filled (hollow) symbols represent p xx (p xy ). 



ers ^(J. More importantly, Fig. 3(b) further reveal that 
the counterflow p xy at v = 1 remains vanishingly small 
at this T for all values of Ap. The data of Fig. 3(a,b) 
illustrate that the dissipation in counterflow transport at 
v = 1 dramatically changes when the layers are imbal- 
anced while the particle- vacancy pairing remains strong, 
as evinced by the vanishing counterflow p xy at v = 1. In 
Fig. 3(c) we summarize our counterflow p xx and p xy vs 
T -1 data, taken at different values of Ap and at constant 
p = 5.45 x 10 10 cm~ 2 . The data substantiate the findings 
of Fig. 3(a,b), namely that the counterflow p xx is very 
sensitive to layer density imbalance, while the counter- 
flow p xy is not. 

Why does a rather small change in layer density dis- 
tribution so dramatically change the counterflow p xx at 
v = 1, while p xy remains vanishingly small? ft is unlikely 
that this behavior results from the small changes in the 
intra-layer interaction energies of the two layers. Disor- 
der, on the other hand, can be the culprit. When Ap < 0, 
the (bottom) layer, which we probe, has a smaller den- 
sity than in the balanced case and is thus more prone to 
disorder since, e.g., the screening of the ionized impurity 



potential is less effective. It is therefore not surprising 
that its p xx at v = 1 is larger than in the balanced case. 
The converse is true when Ap > 0. The data of Fig. 3 
also provide a natural clue for understanding the results 
of Fig. 1. Lowering the density and therefore cL/Ib in the 
balanced case has two consequences which influence p xx 
in opposite directions. First, it strengthens the inter- 
layer interaction and leads to a lowering of p xy , as ob- 
served. Theoretically, we would also expect a reduction 
in p xx with decreasing d/ls- But a second consequence of 
lowering the density is to enhance the effective disorder 
in the bilayer system and thus increase p xx . We believe 
these two, compensating effects are responsible for p xx 
in Fig. 1 being essentially independent of d/ls while p xy 
significantly decreases with decreasing d/ls- The com- 
bination of the data of Figs. 1 and 3 therefore suggests 
that the sample disorder, and the ensuing mobile vortices 
H> 113 1 are the likely culprits for the finite counterflow 
dissipation in our samples. 
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